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ABSTRACT 


Exoplanetary science is increasingly prioritizing efforts toward direct imaging of planetary systems, 
with emphasis on those that may enable the detection and characterization of potentially habitable ex- 
oplanets. The recent 2020 Astronomy and Astrophysics decadal survey recommended the development 
of a space-based direct imaging mission that has subsequently been referred to as the Habitable Worlds 
Observatory (HWO). A fundamental challenge in the preparatory work for the HWO search for exo- 
Earths is the selection of suitable stellar targets. Much of the prior efforts regarding the HWO targets 
has occurred within the context of exoplanet surveys that have characterized the stellar properties for 
the nearest stars. The preliminary input catalog for HWO consists of 164 stars, of which 30 are known 
exoplanet hosts to 70 planets. Here, we provide a dynamical analysis for these 30 systems, injecting a 
terrestrial planet mass into the Habitable Zone (HZ) and determining the constraints on stable orbit 
locations due to the influence of the known planets. For each system, we calculate the percentage 
of the HZ that is dynamically viable for the potential presence of a terrestrial planet, providing an 
additional metric for inclusion of the stars within the HWO target list. Our analysis shows that, for 11 
of the systems, less than 5096 of the HZ is dynamically viable, primarily due to the presence of giant 
planets whose orbits pass near or through the HZ. These results demonstrate the impact that known 
system architectures can have on direct imaging target selection and overall system habitability. 


Keywords: astrobiology — planetary systems — planets and satellites: dynamical evolution and stability 


1. INTRODUCTION 


Although thousands of exoplanets have now been dis- 
covered, a full census and exploration of nearby plan- 
etary systems has yet to be realized. One of the pri- 
mary reasons for incompleteness regarding nearby plan- 
etary architectures are the observational biases inherent 
within the two detection techniques responsible for the 
vast majority of current exoplanet discoveries: the ra- 
dial velocity (RV) and transit methods (Cumming et al. 
2008; Kane & von Braun 2008; Ford 2014; Winn & Fab- 
rycky 2015; Kipping & Sandford 2016; He et al. 2019). 
Both of these methods are heavily biased toward large, 
close-in planets, with strong dependencies on survey de- 
sign, data precision, observing strategy, and the time 
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baseline of observations (Kane et al. 2007; Ford 2008; 
von Braun et al. 2009; Wittenmyer et al. 2013; Barclay 
et al. 2018). However, RV surveys that have been op- 
erating for multiple decades have been able to probe 
beyond the snowline, providing insight into the preva- 
lence of giant planets, even at relatively large separa- 
tions from the host star (Wittenmyer et al. 2020; Fulton 
et al. 2021; Rosenthal et al. 2021). The architectures of 
these systems are informative for their orbital dynamics 
and evolution, that can also be used to constrain the 
possible presence of additional planets within those sys- 
tems (Barnes & Raymond 2004; Kopparapu & Barnes 
2010; Kane 2015; Kane & Blunt 2019; Kane 2023; Kane 
& Fetherolf 2023). Of particular interest is the potential 
for harboring planets within the Habitable Zone (HZ) of 
the system (Kasting et al. 1993; Kane & Gelino 2012; 
Kopparapu et al. 2013; Kane 2014; Kopparapu et al. 
2014; Kane et al. 2016; Hill et al. 2018; Kane et al. 2020a; 
Hill et al. 2023). Such HZ planets are expected to be 
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the target of extensive follow-up observations, especially 
those nearby systems whose host star brightness and 
angular separation of their planets may enable direct 
imaging as a means to characterize the planetary atmo- 
spheres (Brown 2015; Barclay et al. 2017; Kane et al. 
2018; Kopparapu et al. 2018; Stark et al. 2020; Li et al. 
2021). 

Among the recommendations of the 2020 Astronomy 
and Astrophysics decadal survey (hereafter Astro2020) 
was the prioritization of a space-based mission with the 
capability to directly image terrestrial planets within the 
HZ !. The currently adopted name for this recommended 
mission is the Habitable Worlds Observatory (HWO), 
and considerable efforts are being undertaken to define 
the mission designs that would maximize the science 
yield for terrestrial planet characterization (Vaughan 
et al. 2023; Stark et al. 2024). Work is also being carried 
out to both define and characterize the target stars that 
are best suited for HWO observations. An initial list of 
164 HWO target stars has been provided to the broader 
community for assessment and discussion (Mamajek & 
Stapelfeldt 2024), the stellar properties of which have 
been further refined by Harada et al. (2024) and ex- 
amined in the context of yield simulations by Tuchow 
et al. (2024). Of these 164 nearby stars, all of which 
lie within 24 pcs, 30 are currently known to host exo- 
planets, with a total inventory of 70 known exoplanets. 
These exoplanet detections have primarily been enabled 
through long-term RV observations, revealing the Kep- 
lerian orbits and architectures down to the RV precision 
of the corresponding surveys. Though the architectures 
for these planetary systems are likely incomplete, even 
partial knowledge of the planets present enables a dy- 
namical analysis that can inform the likelihood that each 
system has a HZ that contains dynamically viable or- 
bital locations for terrestrial planets. Given the need to 
determine a high quality list of HWO targets, depriori- 
tizing those systems that are dynamically unsuitable is 
a crucial step in the target assessment process. 

In this work, we present the results of a dynamical 
study for the known exoplanet hosts within the cur- 
rent proposed HWO target list provided by Mamajek 
& Stapelfeldt (2024), that specifically provides an as- 
sessment for the dynamical stability within the HZ of 
those systems. Section 2 provides a description of the 
HWO stellar sample, with a particular emphasis on the 
known exoplanetary systems within that list. The de- 
scription includes a tabulation of the stellar and plan- 
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etary properties for the known exoplanet systems, and 
the methodology for the calculation of the system HZ 
boundaries. The structure and results from the dynam- 
ical simulations are described in Section 3, with a quan- 
tification of the available HZ that is dynamically stable. 
We discuss the implications and limitations of our work 
in Section 4, and provide concluding remarks and sug- 
gestions for future work in Section 5. 


2. THE HABITABLE WORLDS OBSERVATORY 
SAMPLE 


Here, we provide the details of the potential HWO 
target list that harbor known exoplanets, and discuss 
the calculation of the HZ for these systems. 


2.1. Known Planetary Systems 


From the 164 stars identified by Mamajek & 
Stapelfeldt (2024), we extracted the 30 stars that are 
currently known to host exoplanets (Akeson et al. 2013). 
'The full list of planets for all of the included systems, 
including their relevant stellar and planetary properties, 
are provided in Table 1. The selection of relevant pa- 
rameters means those that are needed for the calculation 
of the HZ and/or the dynamical simulations, described 
in Section 3. Note that the HH Peg (HD 206860) system 
was not included in our analysis since the companion is 
a brown dwarf at very wide separation from the host 
star, detected using the direct imaging technique (Luh- 
man et al. 2007). The sources for the stellar/planetary 
parameters used in all cases are included in the final 
column of Table 1, where the source selection is based 
upon the most complete data source. 

'The stellar multiplicity of a planetary system can have 
a significant dynamical effect on HZ planets (Holman 
& Wiegert 1999; Kane & Hinkel 2013; Simonetti et al. 
2020). Of the 30 known exoplanet hosts in our sample, 
9 of the stars are known to have stellar or brown dwarf 
companions. The HD 3651 system has a brown dwarf 
with a projected separation of 476 AU away the primary 
(Luhman et al. 2007). The HD 9826 system contains a 
red dwarf star with a projected separation of 750 AU 
from the primary (Mason et al. 2001). HD 26965 is a 
triple system, whereby the B and C components are a 
white dwarf and red dwarf star, respectively, and or- 
bit each other at a distance of ~400 AU from the pri- 
mary (Bond et al. 2017). The HD 75732 system con- 
tains a red dwarf star with a projected separation of 
1100 AU from the primary (Marcy et al. 2002). The 
HD 102365 includes a red dwarf star with a projected 
separation of 211 AU from the primary (Raghavan et al. 
2010). HD 1154044 is in a binary orbit with a red dwarf 
star that has projected separation of 289 AU (Alonso- 
Floriano et al. 2015). HD 147513 is in a binary orbit 
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with a white dwarf separated by ~5300 AU (Porto de 
Mello & da Silva 1997). The HD 190360 system includes 
a red dwarf star that is separated by at least 1266 AU 
from the primary (Naef et al. 2003) and is probably at a 
much larger distance (Feng et al. 2021). The HD 209100 
system contains a brown dwarf binary pair that is sepa- 
rated from the primary by 1500 AU (Scholz et al. 2003; 
Luhman et al. 2007; King et al. 2010). Given the wide 
separations for all of these binary/multiple systems, and 
the considerable uncertainties in their orbital parame- 
ters, the stellar or brown dwarf companions were not 
included in our subsequent dynamical analysis. 

There are several further important notes regarding 
the included systems. The vast majority of the 70 plan- 
ets are not known to transit their host star, and so 
the planetary masses used are minimum masses in most 
cases, the implications of which are discussed in Sec- 
tion 4. Exceptions to this include HD 39091 c (Huang 
et al. 2018), HD 75732 e (Demory et al. 2011; Kane 
et al. 2011; von Braun et al. 2011; Winn et al. 2011), 
HD 136352 b, c, and d (Kane et al. 2020b; Delrez et al. 
2021), and HD 219134 b and c (Gillon et al. 2017), all of 
which are known to transit their host stars. There are 
also planets in our list for which orbital inclination mea- 
surements have been made. For HD 39091, planet b has 
an inclination of ~50°, resulting in a mutual inclination 
of ~40° with respect to planet c (Xuan & Wyatt 2020), 
and we assume planet d is coplanar with planet b. For 
HD 115404A, planet c has an inclination of ~25°, and 
we assume planet b is coplanar (Feng et al. 2022). For 
HD 140901, planet c has an inclination of ~170°, and 
we assume planet b is coplanar Feng et al. (2022). For 
HD 160691, astrometry places an upper limit on the in- 
clination of ~60° for the known planets (Benedict et al. 
2022). For HD 209100 b, Feng et al. (2019) provides an 
inclination of ~64°, thus establishing the actual plane- 
tary mass. Planet masses reported in the literature as 
Earth mass units (Ma) have been converted to Jupiter 
mass units (Mz) in Table 1 for consistency. For cases 
where the eccentricity is reported as zero, the argument 
of periastron is set to 90°, corresponding to the location 
of inferior conjunction. Finally, note that the existence 
of the planet orbiting HD 26965 (Ma et al. 2018) has 
been disputed from RV follow-up data (Burrows et al. 
2024), but we elected to err on the side of caution and 
included the planet in our list. 


'The exoplanet population included in our sample is 
represented in Figure 1, showing the distributions of 
semi-major axes and eccentricities. The ranges of these 
values within the sample are 0.015-11.6 AU and 0.0- 
0.64 for the semi-major axes and eccentricities, respec- 
tivelv. The data point sizes are loearithmicallv propor- 
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Figure 1. Distribution of semi-major axes and eccentricities 
for all planets included in our sample. The size of the plotted 
data are logarithmically proportional to the planet mass. 


tional to the planet mass, with a range of 0.0055 M; 
(tau Ceti g) to 12.6 M; (pi Men b). There are several 
features to consider regarding the parameter distribu- 
tion shown in Figure 1. There is a clear trend toward 
higher mass planets at larger separations. For exam- 
ple, the mean planet mass for semi-major axes less than 
1 AU is 0.19 M;, whereas the mean planet mass beyond 
1 AU is 1.19 M;. A significant factor in this planetary 
mass trend is the observational bias inherent in the RV 
method, whereby the RV amplitude decreases with in- 
creasing semi-major axis (Cumming 2004; O’Toole et al. 
2009). A similar bias impacts the observed eccentricity 
distribution, since such planets produce higher RV am- 
plitudes, albeit during a limited period of the orbital 
phase (Kane et al. 2007; Shen & Turner 2008; Zakam- 
ska et al. 2011; Kane et al. 2012). Moreover, for those 
long-period planets with insufficient phase coverage, the 
relatively slow apastron passage may result in eccentric- 
ity misclassifications (Anglada-Escudé et al. 2010; Hara 
et al. 2019; Wittenmyer et al. 20192). Even so, the com- 
bination of relatively high masses and eccentricities pro- 
duces numerous planetary architectures that exhibit sig- 
nificant dynamical consequences within the HZ of their 
systems. 
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Table 1. Stellar and planetary parameters for HWO known systems. 


Star NP! M, Ter Planet Mj a e w Source 
(Mo) (K) (Mz) (AU) (°) 
HD 3651 (54 Piscium) 1 0.799 5221 b 0.228 0.295 0.645 243  Wittenmyer et al. (2019b) 
HD 9826 (ups And) 3 1.29 6156 b 0.675 0.05914 0.0069 0 Rosenthal et al. (2021) 
c 1.965 0.8265 0.266 58.2 
d 4.1 2.517 0.294 73.8 
HD 10647 (q1 Eridani) 1 1.11 60218 b 0.94 2.015 0.15 212 Marmier et al. (2013) 
HD 10700 (tau Ceti) 4 0.78 5333 g 0.0055 0.133 0.06 395.34 Feng et al. (2017b) 
h 0.0058 0.243 0.23 7.45 
e 0.0124 0.538 0.18 22.35 
f 0.0124 1.334 0.16 119.75 
HD 17051 (iota Hor) 1 1.34 6167 b 2.27 0.92 0.14 309 Stassun et al. (2017) 
HD 20794 (82 Eridani) 4 0.813 5401 b 0.0089 0.127 0.27 386.17 Feng et al. (2017a) 
c 0.0079 0.225 0.17 31.51 
d 0.0111 0.364 0.25 252.67 
e 0.015 0.509 0.29 268.14 
HD 22049 (eps Eri) 1 0.81 5020 b 0.651 3.5 0.044 350 Rosenthal et al. (2021) 
HD 26965 (40 Eri) 1 0.78 5072 b 0.0267 0.219 0.04 148.969 Ma et al. (2018) 
HD 33564 (kappa Cam) 1 1.25 6250 b 9.1 1.1 0.34 205 Galland et al. (2005) 
HD 39091 (pi Men) 3 1.07 5998 c 0.0114 0.06805 0 90 Hatzes et al. (2022) 
d 0.0421 0.499 0.22 323 
b 12.6 3.2826 0.6396 331.03 
HD 69830 3 0.86 5385 b 0.0321 0.0785 0.1 340  Lovis et al. (2006) 
c 0.0371 0.186 0.13 221 
d 0.057 0.63 0.07 224 
HD 75732 (55 Cancri) 5 0.905 5172 e 0.0251 0.01544 0.05 86 Bourrier et al. (2018) 
b 0.8036 0.1134 0 90 
c 0.1611 0.2373 0.03 2.4 
f 0.1503 0.7708 0.08 262.4 
d 3.12 5.957 0.13 290.9 
HD 95128 (47 UMa) 3 1.06 5872 b 2.53 2.1 0.032 334 Gregory & Fischer (2010) 
c 0.54 3.6 0.098 295 
d 1.64 11.6 0.16 110 
HD 95735 2 0.3899 3712 b 0.0085 0.07879 0.063 330 Hurt et al. (2022) 
c 0.0428 2.94 0.132 63 
HD 102365 1 0.85 5630 b 0.0503 0.46 0.34 105 Tinney et al. (2011) 
HD 114613 1 1.27 5641 b 0.357 5.34 0.458 196 Luhn et al. (2019) 
HD 115404A 2 0.83 5019 b 0.097 0.088 0.232 259 Feng et al. (2022) 
10.319 11.364 0.211 312 
HD 115617 (61 Vir) 3 0.942 5577 b 0.016 0.050201 0.12 105 Vogt et al. (2010) 
c 0.057 0.2175 0.14 341 


Table 1 continued 
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Table 1 (continued) 


Star NPi M, Teg Planet M3 a e w Source 
(Mo) (K) (Ms) | (AU) (°) 
d 0.072 0.476 0.35 314 
HD 136352 (nu2 Lupi) 3 0.87 5664 b 0.0149 0.0964 0 90 Delrez et al. (2021) 
c 0.0354 0.1721 0 90 
d 0.0278 0.425 0 90 
HD 141004 (lambda Ser) 1 1.05 5885 b 0.0428 0.1238 0.16 19.48 Rosenthal et al. (2021) 
HD 140901 2 0.99 5586 b 0.0503 0.085 0.472 108 Feng et al. (2022) 
c 6.284 7.421 0.607 298 
HD 143761 (rho CrB) 4 0.95 5817 e 0.0119 0.1061 0.126 359.4 Brewer et al. (2023) 
b 1.093 0.2245 0.038 269.64 
c 0.0887 0.4206 0.096 9.7 
d 0.068 0.827 0 90 
HD 147513 (62 G. Sco) 1 1.11 5883 b 1.21 1.32 0.26 282 Mayor et al. (2004) 
HD 160691 (mu Ara) 4 1.13 5773 d 0.033 0.0923 0.16 197 Benedict et al. (2022) 
e 0.439 0.9296 0.091 193 
b 1.665 1.5224 0.036 39 
1.873 5.0937 0.022 84 
HD 189567 2 0.83 5726 b 0.0267 0.111 0 90 Unger et al. (2021) 
c 0.022 0.197 0.16 219 
HD 190360 2 0.99 5537 c 0.0677 0.1294 0.165 322 Rosenthal et al. (2021) 
b 1.492 3.955 0.3274 13.9 
HD 192310 2 0.8 5166 b 0.0532 0.32 0.13 173 Pepe et al. (2011) 
c 0.076 1.18 0.32 110 
HD 209100 (eps Indi A) 1 0.754 4611 b 3.25 11.55 0.26 77.83 Feng et al. (2019) 
HD 217987 2 0.489 3688 b 0.0132 0.068 0 90 Jeffers et al. (2020) 
c 0.0239 0.12 0 90 
HD 219134 6 0.794 4913 b 0.012 0.038474 0 90 Vogt et al. (2015) 
c 0.011 0.064816 0 90 
f 0.028 0.14574 0 90 
d 0.067 0.23508 0 90 
g 0.0346 0.3753 0 90 
h 0.34 3.11 0.06 215 


T Number of planets. 


Minimum mass, except for 55 Cancri e and the HD 136352 planets. 


boundaries of the HZ are derived from 1D Earth-based 
climate models that calculate the radiative balance for 
which surface liquid water is retained. This approach di- 
2.2. System Habitable Zones vides the HZ into two primary regions: the conservative 
HZ (CHZ) and the optimistic HZ (OHZ). The CHZ is 
bounded by the runaway greenhouse criteria at the inner 
edge, and the maximum CO» greenhouse criteria at the 
outer edge (Kane et al. 2016). The OHZ is an empirical 


'The extent of the HZ for each system was calculated 
using the stellar parameters from the sources provided 
in Table 1. Specifically, we adopt the methodology pro- 
vided by Kopparapu et al. (2013, 2014), in which the 
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Figure 2. HZ and planetary orbits for four of the systems in our sample: HD 10647 (top-left), HD 39091 (pi Men; top-right), 
HD 75732 (55 Cancri; bottom-left), HD 95128 (47 Uma; bottom-right). The orbits are labeled by planet designation. The 
extent of the HZ is shown in green, where light green and dark green indicate the CHZ and OHZ, respectively. 


extension to the CHZ based upon assumptions regard- 
ing retention of surface liquid water in the Venusian and 
Martian evolutionary histories (Baker 2001; Kane et al. 
2014; Way et al. 2016; Orosei et al. 2018; Kane et al. 
2019; Kane & Byrne 2024). The CHZ and OHZ bound- 
aries for each system are provided in Table 2. 

Figure 2 and Figure 3 show top-down views for 8 of 
the planetary systems, that demonstrate the diversity of 
architectures in our known exoplanet host sample. The 
extent of the CHZ and OHZ are shown in light green 
and dark green, respectively. In most of these cases, the 


planetary orbit that is present within the HZ is a gi- 
ant planet with significant perturbation potential. The 
cases where eccentric orbits lie within and/or cross the 
HZ, such as HD 39091 (pi Men), HD 114613, HD 147513, 
and HD 192310, are particularly devastating to orbital 
stability in that region due to angular momentum trans- 
fer that excites other planetary orbits into high eccen- 
tricity regimes (Kane & Raymond 2014). The others 
cases shown exhibit planetary orbits that are more cir- 
cular, and can often allow regions of orbital stability 
at the edges of the HZ. These scenarios may be tested 
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Figure 3. HZ and planetary orbits for four of the systems in our sample: HD 114613 (top-left), HD 147513 (top-right), 
HD 160691 (bottom-left), HD 192310 (bottom-right). The orbits are labeled by planet designation. The extent of the HZ is 
shown in green, where light green and dark green indicate the CHZ and OHZ, respectively. 


with dynamical simulations that explore particle injec- 
tion throughout the HZ region. 


3. ORBITAL DYNAMICS 
3.1. Simulation Design 


We conducted an exhaustive series of dynamical simu- 
lations to test the viability of terrestrial planets orbiting 
within the HZ of the target systems. The simulations 
were conducted using the REBOUND N-body integra- 
tor package (Rein & Liu 2012) that applies the sym- 
plectic integrator WHFast (Rein & Tamayo 2015). Our 


simulations explored the effect and orbital evolution of 
an additional Earth-mass planet in a circular orbit that 
is coplanar with the other planets within each system, 
where 10 equally spaced mean anomalies were chosen for 
the injected planet. For the case of HD 39091 (pi Men), 
where there is a substantial mutual inclination between 
the planets, we injected the additional planet in an or- 
bit coplanar with planet b, since that is the planet that 
dominates the gravitational influence on the HZ for that 
system (see Section 2.1 and Figure 2). The semi-major 
axes explored spanned 20% less than the inner edge 
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of the OHZ to 20% larger than the outer edge of the 
OHZ (see Section 2.2). This region within each system 
was divided into 1000 equally spaced semi-major axis 
starting locations for the injected planet to ensure suf- 
ficiently dynamical sampling of the HZ. Each of these 
semi-major axis and mean anomaly starting locations 
were integrated for 10" years, resulting in 10,000 simula- 
tions performed for each system. The time step for each 
integration was selected based upon the orbital period 
for the innermost body within the system divided by 20, 
where the innermost body is either a known planet or 
the injected planet. 

At each of the simulation starting locations, our sim- 
ulations provided the percentage of the total simulation 
duration for which the planet was able to survive, where 
non-survival means that the injected planet was cap- 
tured by the gravitational well of the host star or ejected 
from the system. From these results, we calculated the 
dynamically viable HZ (DVHZ), which is the percent- 
age of the OHZ locations for which the injected planet 
survived the full 107 year simulation. 


3.2. Habitable Zone Stability Results 


'The methodology described in Section 3.1 was applied 
to each of the systems shown in Table 1. The time re- 
quired to complete the suite of simulations for each sys- 
tem varied enormously, depending upon the number of 
planets in the system, their locations, and the extent 
of the instability induced within the HZ. For example, 
the 55 Cancri system (bottom-left system depicted in 
Figure 2) contains 5 known planets, including one with 
a 0.74 day orbit, resulting in the need for high time 
resolution within the dynamical simulation. By con- 
trast, systems with a single, eccentric giant planet with 
high scattering potential (such as the HD 147513 sys- 
tem) concluded the simulation suite relatively quickly. 

'The simulation results for 8 of the systems are repre- 
sented in Figure 4 and Figure 5, corresponding to the 
same systems shown in Figure 2 and Figure 3, respec- 
tively. The plots provide the percentage of the full 107 
year simulation for which the injected terrestrial planet 
was able to maintain a stable orbit for each semi-major 
axis location. As for Figure 2 and Figure 3, the CHZ 
and OHZ are shown as light green and dark green shaded 
regions, respectively. In all of the 8 cases shown, there 
are significant areas of instability within the HZ result- 
ing from the known planets in the system. Most no- 
table is the case of pi Men (second panel; Figure 4) in 
which the entire region explored is unstable due to the 
incursion of the highly eccentric b planet into the HZ. 
The HD 147513 system (second panel; Figure 5) is sim- 
ilarly affected, although stable regions exist just outside 


of the HZ. The remaining examples shown have varying 
amounts of stable locations within the HZ, including re- 
gions of stability/instability resulting from a combina- 
tion of mean motion resonance (MMR) effects and the 
numerical randomness intrinsic to the simulations (Ray- 
mond et al. 2008; Petrovich et al. 2013; Hadden 2019). 
For example, the results for HD 10647 (first panel; Fig- 
ure 4) exhibits numerous stability islands induced by the 
single giant planet in the system. The 3:2 MMR loca- 
tion at 1.53 AU creates an island of stability whereby 
a planet located there has the potential to avoid grav- 
itational perturbing interactions with the known giant 
planet, similar to the 2:3 MMR of Neptune with Pluto 
(Williams & Benson 1971). Of particular note is the 
5:2 MMR location at 1.09 AU that produces a region 
of instability that coincides with the inner boundaries 
of the OHZ and CHZ. Multi-planet systems can exhibit 
significantly complicated stability patterns due to the 
multitude of MMR locations cascading throughout the 
system, as shown in the case of HD 75732 (55 Cancri; 
third panel; Figure 4). 


To quantify the dynamical simulation results on a 
system-by-system basis, we calculate and provide the 
DVHZ in each case, as described in Section 3.1. These 
calculations are shown in Table 2, along with the inner 
and outer boundaries for the CHZ and OHZ, described 
in Section 2.2. For 11 of the systems, DVHZ is less than 
5096, and there are 4 cases where the DVHZ is 096, in- 
cluding the previously discussed systems of HD 39091 
(pi Men) and HD 147513. In general, the architectures 
for the 11 cases where DVHZ is less than 5096 are domi- 
nated by a giant planet in an eccentric orbit that passes 
near or through the HZ. 

On the other hand, there are 9 cases for which the 
DVHZ is 10096. Cross-referencing these systems with 
Table 1 reveals that their architectures are dominated 
by low-mass, short-period planets with limited gravita- 
tional influence within the HZ. Shown in Figure 6 is 
a histogram of the DVHZ values that summarizes the 
findings of our dynamical simulations. This shows that 
half of the studied systems have DVHZ values that are 
greater than 8096, ensuring that there is substantial dy- 
namically viable space within the HZ of those systems 
based on the currently known architectures, and there- 
fore suitable for potential direct imaging follow-up ob- 
servations. 


4. DISCUSSION 


Given the number of known exoplanetary systems, in- 
cluding those resulting from numerous decades of RV 
observations around bright stars, it is little surprise that 
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Figure 4. Percentage of the simulation that the injected Earth-mass planet survived as a function of the initial semi-major 
axis, where the CHZ is shown in light green and the OHZ is shown in dark green. Shown, from top to bottom, are the results 
for HD 10647, HD 39091 (pi Men), HD 75732 (55 Cancri), and HD 95128. 
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Figure 5. Percentage of the simulation that the injected Earth-mass planet survived as a function of the initial semi-major 
axis, where the CHZ is shown in light green and the OHZ is shown in dark green. Shown, from top to bottom, are the results 
for HD 114613, HD 147513, HD 160691, and HD 192310. 
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Table 2. HZ boundaries and DVHZ results. 


Star OHZin CHZin CHZou OHZour DVHZ 
(AU) (AU) (AU) (AU) (7) 

HD 3651 0.56 0.71 1.27 1.34 92.41 
HD 9826 1.35 1.71 3.00 3.16 0.00 
HD 10647 0.87 1.10 1.92 2.02 38.65 
HD 10700 0.55 0.69 1.24 1.31 . 20.75 
HD 17051 0.94 1.20 2.09 2.20 55.02 
HD 20794 0.72 0.91 1.63 1.72 0.00 
HD 22049 0.45 0.57 1.03 1.09 100.00 
HD 26965 0.52 0.66 1.20 1.26 100.00 
HD 33564 1.58 2.00 3.49 3.68 61.50 
HD 39091 0.95 1.20 2.11 2.23 0.00 
HD 69830 0.60 0.75 1.35 1.42 84.48 
HD 75732 0.59 0.74 1.33 1.41 56.99 
HD 95128 0.94 1.19 2.09 2.20 27.59 
HD 95735 0.13 0.16 0.30 0.32 100.00 
HD 102365 0.57 0.72 1.28 1.85 78.21 
HD 114613 1.54 1.96 3.46 3.65 33.43 
HD 115404A 0.43 0.54 0.99 1.04 100.00 
HD 115617 0.68 0.86 1.53 1.61 44.93 
HD 136352 0.77 0.97 1.72 1.82 100.00 
HD 140901 0.69 0.87 1.54 1.62 89.06 
HD 141004 1.08 1.36 2.40 2.53 100.00 
HD 143761 1.01 1.28 2.25 2.38 | 88.17 
HD 147513 0.74 0.93 1.64 1.73 0.00 
HD 160691 1.00 1.26 2.23 2.35 5.39 
HD 189567 1.09 1.38 2.45 2.58 100.00 
HD 190360 0.82 1.04 1.84 1.04 85.67 
HD 192310 0.48 0.61 1.10 116 19.14 
HD 209100 0.38 0.48 0.87 0.92 100.00 
HD 217987 0.16 0.20 0.38 0.40 100.00 
HD 219134 0.41 0.52 0.95 1.00 99.26 


the HWO target list contains known exoplanet systems. 
The description of these 30 systems in Section 2 shows 
that the architectures of these systems is diverse, in- 
cluding low-mass multi-planet systems and high-mass 
single-planets in eccentric orbits. The latter category in- 
cludes giant planets beyond the snow line, the frequency 
of which appears to be relatively rare (Albrow et al. 
2001; Wittenmyer et al. 2011, 2016; Fulton et al. 2021; 
Rosenthal et al. 2021; Bonomo et al. 2023), and can 
greatly influence the volatile delivery within the system 
(Raymond & Bonsor 2014; Raymond & Izidoro 2017; 
Venturini et al. 2020; Kane & Wittenmyer 2024). The 
high eccentricity cases may be the result of planet-planet 
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Figure 6. Histogram of the DVHZ calculations for all 30 of 
the known systems studied in this work. 


scattering events (Carrera et al. 2019) and can dynam- 
ically eliminate the possibility of further planets being 
present in the system (Brewer et al. 2020). The grow- 
ing understanding of such dynamically volatile systems, 
including their formation and prevalence, will greatly 
inform our ability to infer the possible presence of HZ 
exoplanets (Barnes & Raymond 2004). 

As described in Section 3.1, the integration time for 
each of our simulations was 10" years. The integration 
time was chosen to minimize the effects of orbital ele- 
ment error propagation (Ford 2005; Kane et al. 2009), 
and to create a feasible pathway to explore the full pa- 
rameter space for all 30 targets, each of which required 
thousands of simulations. However, chaotic systems are 
often the result of divergent orbital eccentricities that 
can require significant time to develop. It is thus pos- 
sible that, even for locations identified as stable in our 
simulations, instability occurs beyond the 10" year in- 
tegration time (Gozdziewski et al. 2001). Although our 
criteria for stability is a first-order approach that neces- 
sitates all planets survive the duration of the dynamical 
simulation, this criteria has successfully been validated 
against chaos indicators (Dvorak et al. 2010) and applied 
to numerous exoplanetary systems (Menou & Tabachnik 
2003; Dvorak et al. 2003; Kane 2016). 

As noted in Section 2.1, the preponderance of plan- 
ets in our sample are not known to transit their host 
stars, whereby the true mass of the planet may be re- 
alized. Thus, the planetary masses used are generally 
minimum masses. A consequence of this knowledge lim- 
itation is that the dynamical influence of the planets on 
each other within a system, and their effects on the in- 
jected planet, may be larger than that calculated and 
presented here. For that reason, the results presented 
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in Section 3.2 can be considered a lower limit on the 
HZ orbital instability. Furthermore, the observational 
bias in the RV and transit methods mentioned in Sec- 
tion 2.1 means that the architectures of the known sys- 
tems are likely incomplete, with additional planets yet 
to be discovered in these systems. However, it is worth 
repeating here that we assumed coplanar orbits for the 
injected planets in our simulations (coplanar with planet 
b in the case of pi Men). A full exploration of terrestrial 
HZ planets whose orbits are mutually inclined from the 
known planets may yet reveal stable locations that are 
otherwise unviable. 


5. CONCLUSIONS 


The effective transition of the exoplanet community to 
prioritizing direct imaging techniques requires that the 
data for the expected stellar targets be sufficiently lever- 
aged toward a robust target list. Of the 164 targets pro- 
vided by Mamajek & Stapelfeldt (2024) for direct imag- 
ing with HWO, only ~20% are presently known to har- 
bor exoplanets. Most, if not all, of the remaining stars 
likely also have planetary companions, but in most cases 
the current database of RV observations may be suffi- 
cient to rule out the presence of giant planets in or near 
the HZ. For those that are known exoplanet hosts, our 
results show that 11 of the systems have less than 50% 
of the HZ that is dynamically viable (DVHZ < 50%), 
including 4 cases where there is no dynamically possible 
orbit (ups And, eps Eri, pi Men, and 62 G. Sco). These 
systems should therefore be treated with extreme cau- 
tion when considered as possible HWO targets. How- 
ever, over half of the tested systems have DVHZ values 
greater than 80%, making them suitable HWO targets 
from the perspective of orbital dynamics. Moreover, our 
overall results are not necessarily representative of what 
may be expected for the broader list of HWO targets 
where no exoplanets have yet been detected. Giant plan- 
ets are relatively rare among the exoplanet population, 
and the detection methods utilized are biased toward de- 
tecting those kinds of planets more quickly than lower- 
mass planetary counterparts. This means that current 
exoplanet surveys may have already detected many of 


the nearby systems that contain significant perturbing 
agents (giant planets) for possible terrestrial HZ planets. 

With respect to exoplanet architectures, the refine- 
ment of the HWO target sample will greatly benefit 
from the continued search for exoplanets. The era of 
extreme precision RVs (Fischer et al. 2016) provides op- 
portunities to detect previously undiscovered exoplan- 
ets around the nearby stars. These exoplanet detections 
will include low-mass terrestrial exoplanets and high- 
mass planets with low-inclinations relative to the plane 
of the sky. Additional planets are also likely to be found 
in the known systems, further constraining the stability 
profiles for the HZ region, as described in this work. Re- 
sults from the Gaia mission (Gaia Collaboration et al. 
2021) will aid in the detection of long-period giant plan- 
ets and allow the true mass of known RV giant planets to 
be determined. All of these precursor observational ef- 
forts will continue to form an essential component for the 
success of the Astro2020 recommendations, and maxi- 
mize the yield of HWO and the characterization of HZ 
terrestrial planets. 
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